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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2768 


SUPERSONIC FLOW WITH WHIRL AND VORTICITY IN AXISYMMETRIC 

CHANNELS 

By Ralph <1. Eschborn 


SUMMARY . . . 

The axially symmetric supersonic steady flow is treated for a non- 
viscous fluid flowing with whirl component between two arbitrary coaxial 
surfaces of revolution. The equations that describe this motion Eire 
expressed in characteristic coordinates and in this form are used to 
determine the meridional (axial- radial) velocities for arbitrary distri- 
butions of the tangential velocity components. 

Solution of the supersonic characteristic equations is investigated 
for arbitrary axisymmetric flow fields with and without vorticity. This 
solution includes cases in which the channel shape and the inlet velocity 
distribution are prescribed. Also included are cases in which one of 
the channel surfaces, the desired velocities on that surface (within 
certain limitations), and the inlet velocity distribution are prescribed. 


INTRODUCTION 

A survey of the literature showed that a method of design or analy- 
sis for supersonic flow with whirl in the three-dimensional region of an 
axially symmetric channel formed by two coaxial surfaces of revolution 
without blades is not available. In the past, axially symmetric channels 
of this type have been designed from one-dimensional flow theory or the 
stream-filament method, which does not consider the character of the flow 
(subsonic or supersonic) . More recently, a method of solution was out- 
lined in reference 1 for subsonic flow by means of an iteration process 
using orthogonal curvilinear coordinates for arbitrary known boundaries. 

A method is described herein which was devised at the NACA Lewis 
laboratory for analyzing supersonic flow in this type channel. This 
method can be used to obtain the velocities in the flow field when the 
boundaries are known with either prescribed rotational or irrotational 
flow at the inlet . It can also be used to obtain the opposite boundary 
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to satisfy a prescribed boundary with the desired velocities on that 
boundary (within certain' regions -which are defined later) for either a 
prescribed rotational or irrotational flow at the inlet or the outlet. 

The velocity components of the flow equations are treated directly 
in accordance with the procedure outlined in reference 2. The same 
technique for axisymmetric flow in rotating impellers is utilized in 
reference 3. 

In addition to describing the method in a mathematical outline, 
this report gives numerical examples which demonstrate same of the 
phenomena present in the flow in either supersonic vaneless difusers or 
the channel between the outlet of a supersonic compressor impeller and 
an annular cascade of supersonic diffuser blades. The method is par- 
ticularly suited to determine pressure gradients and possible shock con- 
figurations when supersonic flow exists in radial diffusers and in 
annular channels at the inlet or outlet of a compressor or turbine. 

This allows the determination of possible flow separation within the 
channel and suggests design modifications to eliminate the losses result- 
ing from shocks. 


DERIVATION OF .FLOW EQUATIONS 

The flow equations to be developed consider the steady flow of a 
nonviscous compressible, fluid under conditions of axial symmetry and 
isentropic state changes. The coordinates used and the velocity com- 
ponents chosen are shown in figure 1. 


Stream Function 

With the use of the symbols defined in the appendix, the continuity 
equation for steady flow is 


V • (pV) = 0 

In cylindrical coordinates, for axial symmetry (partial derivations of 
the fluid properties with respect to 9 and zero) , this equation becomes 


(pvr) + (pwr) = 0 


( 1 ) 


♦ 




Equation (l) implies the existence of a stream function (of two variables 
r and z) such that 
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dijr 

= -pvr 

H=° 

In vector notation, equations (2) are 

V+ = prj X V (2a) 

Equation of Motion 

For steady nonviscous flow, the equation of motion is 

(T-V) V - - ® 

For isentropic conditions, the pressure term may he eliminated since 
Hence, 

(V-V) V = - VH 

The general energy equation for the conditions of no work done and no 
heat transfer states that the stagnation or total energy of the fluid 

% = H +\ 

is constant on any streamline. If the gradient is taken, this equation 
becomes 

V 2 

VHrp = VH + V ~ 2 ~ 

Upon substitution into the equation of motion, 

_ _ V 2 

(V-V) V - V -g* = -vap 
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or 


v X (VxV) =7Hr (3) 

which, is one form of the equation of motion} for later development into 
the characteristic equations, however, it iB convenient to transform the 
equation in the following manner: If the scalar product of equation (3) 

with the unit vector j is taken and the fact that for social symmetry 


j • VHj — 0 

is considered, the resulting equation is 


(4) 


V* [j X (VxV)] = 0 (5) 

This equation further reduces to 

V * V (ru) = 0 (6) 

since for axial symmetry, 

j X (VXV) = j V (ru) (7) 

From equations (2a), (3), (4), (6), and (7), it i£ evident that V (ru) , 
VHj, and Vf axe normal to the velocity vector V and the unit vec- 
tor j. Therefore, VHrji and V(ru) must he parallel to V1 r and hence 
Efji and (ru) must he functions of i|r only (reference 4)} that is, 


Cfflrj 

v Hr = vt 


(8) 


v (ru) = v^ 


d(ru) 

di|r 


( 9 ) 


For any vector A and any unit vector n, 
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may "be obtained by expansion of the product n X n X A. Therefore the 
curl VXV can "be written 


VXV = j [j • (VXV)] - i X [j X (7xv)j 


Insert ing this expression for the curl into equation (3) , expanding the 
triple vector product, and utilizing equation (5) give the following 
equation: 


V x j [j • C7XV)] + (J.Y) [i X (7 XV)] = V% 

When equations (2a), (7), (8), and (9) are inserted, this equation 
"becomes 


V* 


fj’faXV) u d(ru) ^t] 

[_ pr r cb{r + dijr J u 


( 10 ) 


which is the most convenient vector form of the equation of notion. 


General Equations in Cylindrical Coordinates 

ft 

The equations for the flow have "been developed, thus far, in vector 
notation. By the use of cylindrical coordinates, the equations are 
translated to scalar notation. 

Continuity equation. - If the vector form of the continuity equa- 
tion (1) is expanded and divided "by p, the following equation is 
obtained: 


4 


V • V + V • = 0 

By use of the isentropic flow relation 

VP VH 
P = a 2 

■the energy equation, and the fact that V-VErp = 0, there results 

- 1 - V 2 

— ^ ^ 2 = 0 

a 
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In cylindrical coordinates, for axial symmetry, this expression becomes 

2 ’ 


^ ' i) Tr + i 1 ' lO 1,2 ' Is ^ + Tz) " (Is) % " (Is 


"s + ? = ° 


From equation (6), the Up and. u z terms may be eliminated, and. the 
continuity equation in cylindrical coordinates becomes 


v2 

1 - — I v„ + 

Q» 


, W 2 \ VW . , ( U 2 \ V , . 

^ - a 2 ) w z “ a 2 ( w r + T z) + i ^- + a 2 ) r 83 ® 


Equation of motion . - From equation (10) , since i|f is not constant, 
the scalar equation of motion becomes 


d(ru) dH I . . 

P u ~W~~ pr W ~ - w r ) = 0 


( 12 ) 


wherein j • fax V ) has been expressed in cylindrical coordinates for 
axial symmetry. Equation (12) is applicable for inlet flow conditions 
which have tangential velocity and nonuniform energy level. If the 
energy level of the flow entering is uniform, the term dHj/dijr vanishes, 

corresponding to the vector equation VX(VXV) = 0, which was solved in 
reference 1 for subsonic flow. 

If in the flow entering (ru) is a constant or zero, the term 
d(ru)/di|f vanishes. If the flow entering is isoenergetic (dH/p/diJf = 0) 
and (ru) is constant or zero, the flow is irrotational and equa- 
tion (12) reduces to the familiar expression for two-dimensional flow, 

bv _ bw 

Sz §r 


Characteristic Equations 

The supersonic flow -under consideration is governed by equa- 
tions (ll) and (12), which form a system of two equations for the two 
variables w and v as functions of both r and z. This system of 
equations can be solved by the method outlined in reference 2 and leads 
to two basic characteristic equations having new independent variables 
t and T), instead of r and z. Each of these two basic characteris- 
tic equations contains derivatives with respect to only one of the 
independent variables £ or rj, instead of both independent variables 
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as in the case of equations (ll) and (12) . The independent variables f 
and T| are the characteristic parameters. Two additional equations are 
also obtained and these are used to find the position coordinates r 
and z} each of these two equations contains derivatives of r and z 
with respect to only one of the independent variables | or q. The 
characteristic lines have slopes dr/dz. The slope of the characteristic 
line designated by £_j_ corresponds to variable | and constant T), and 

the slope designated corresponds to constant | and variable q. 

The results in applying this method are shown in the following equations 
in which £ represents either of the slopes of the characteristic lines 
and a represents either of the coordinates | or T): 

r ff = £ z cr (13) 

and 

* 

W a + Xv a + Yz 0 = 0 (14) 

where 

<15) 

2 

a 



In the terms In which a choice of sign exists in equations (15) to (17)/ 
the upper sign refers to a coefficient with a plus subscript and the 
lower sign to a coefficient with a minus subscript. Thus, X + and Y + 

correspond to variable £ and constant r\ , and X_ and Y_ correspond 

to constant £ and variable tj. Examination of equations (15) to (17 ) 
shows that = X_ and £_ = X + and that the system is hyperbolic if 
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> 1. As shown In reference 2, hodograph solutions tabulated in 

advance of application to the solution of any problem that may arise 
cannot be made, since equation (14) is nonhomogeneous (Y f 0) and 
cannot be integrated directly. As a result, the four equations obtained 
from equations (13) and (14), when a is equal to I or tj, must be 
integrated simultaneously in a step-by-step procedure. 

Although the following equations will not be used in this report, 
they sire included to show a method of solution when the flow is sub- 
sonic. Combining equations (2) and (12) gives' 


*rr + *zz 


_ c) infer), 

5r v r 


5 ln(pr) 

Sz 


♦a + (pr ) 2 


u d(ru) 
_r oi|r 


dErp 

" Sip. 


= o 


( 18 ) 

which is the principal equation for solution of subsonic flow. This 
equation may be shown to be of the elliptic type if the meridional veloc- 
ity V m = Vv2” + w^ is less than the speed of sound. The subsonic flow 

problem may be solved by a system of equations consisting of equa- 
tions (2), (18), and the density relation for isentropic flow 


P 



2 


1 


v 2 y - 1 


(19) 


This system of equations can be solved far given boundary conditions by 
relaxation methods (references 5 and 6) to give the distribution of i|r 
in the meridional plane, from which the distribution of the velocities 
may be obtained from equation (2). The application of relaxation 
methods to a similar system of equations for the subsonic flow through 
radial- and mixed-flow centrifugal compressors is shown in reference 7. 
Another method of solution for a similar system of equations is dis- 
cussed in detail in reference 8 for flow through turbomachines with 
arbitrary hub and casing shapes. — 


NUMERICAL HlOCEDURE 

Two types of problem exist, the direct or analysis problem in which 
the velocities are determined for given boundaries of an arbitrary 
c hanne l shape, and the inverse or design problem in which the boundaries 
are determined for prescribed velocities.’’ Both the supersonic analysis 
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and tbe design problems may be solved by the characteristic equations 
and will be discussed in detail after a consideration of the method of 
computation. In general, suppose that conditions are known at P]_ and 

Pg (fig* 2 ) and that it is desired to deter min e conditions at Pg. 
Equation (13) written in difference form for these points is 


r 3 “ r l = ( Vl-3 (Z 3 - Z l } (20) 

r 3 ” r 2 = (O 2-3 (z 3 “ Z 2 } (21) 

where (^ + )i-3 indicates the average value of between points P^ 

and P and (£ ) 2 3 indicates the average value of £ between 
points P 2 and P 3 . Similarly; equations (14) may be written in dif- 
ference form; 


w 3 ~ V 1 + ( X 4-)l-3 ( v 3 “ T l) + ^+^ 1-3 ( z 3 ~ z l) ~ 0 (22) 

w 3 " v 2 + ( X -)2-3 ( t 3 " v 2^ + CO 2-3 ( z 3 * z 2^ = 0 ( 23 ) 

These equations also indicate that the coefficients are averages of the 
values at the appropriate points. Equations (20) and (21) may be solved 
for r^ and Z 3 ; after which equations (22) and (23) may be solved for 

W 3 and v^. For the first approximation, only the values of the coeffi- 
cients X; Y are known at points P^ and P 2 ; and these values are 

used in place of averages. After the computation of approximate values 
of the position coordinates and the velocity components at P 3 ; the 

coefficients £; X, Y may be computed at this point. Average coeffi- 
cients may now be found and a better approximation to the position 
coordinates and the velocity components obtained at P 3 may be obtained. 

Convergence to an answer of satisfactory accuracy; for a given interval 

or net size, is obtained by repeating this process as often as required. 

The answers will converge to the solution with fewer iterations; how- 
ever ; the smaller the interval between points. The system of equa- 
tions (20) to (23) may also be applied to points p i' p and P^ if 

the known points happen to have their positions related to each other 
in this manner instead of in the manner indicated by P^, P 2 ; and P 3 . 

In order to determine d(ru)/d\Jr and dH^/dilr when they are not equal to 
zero; In each approximation obtained In the preceding process; ^ must 
be computed and can be obtained from 


* 
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+3 " * “ 



(+ r £ + t z ) 


dz = 





dr 


(24) 


In this equation, £ + is used if the integration is performed along the 
constant-^ line, or £_, if along the constant- | line. Since the deriva- 
tives d(ru)/dijr and dBjj/di|f are total derivatives, that is, ru and Bp 

are functions of iff only, these derivatives can he evaluated from the 
inlet flow conditions for any streamline (^-function) . Hence from the 
value of 1^3, either d(ru)/d^, or dELp/ dijr, or hoth may he found at P3. 

Suppose that conditions are known at Pg (fig. 3) and it is desired 
to determine conditions at the intersection P3 of the constant-^ line 
which passes through Pg and the known boundary AB. The boundary AB 
is a known function and may he specified as 


r = g(z) 


In difference, form, equation (13) written for these points along a 
constant- | line is equation (21) and since 

r 3 = 83(2) (25) 

equations (2l) and (25) may he solved for r3 and z 3 . The slope of the 
boundary at point P 3 is given by 


/ dg\ _ / dr\ _ ^3 
\&z/Z \ dz /3 w 3 


(26) 


and since equation (23) for these points along a curve of constant | 
is the difference form of equation (14), equations (23) and (26) may he 
solved for v 3 and w 3 . As previously mentioned, values at the known 

point are used in lieu of averages for the first approximation and the 
process is repeated as often as necessary to obtain the desired accuracy. 
Conditions at point P£ are determined in a similar manner. If either 

d(ru)/d\|r 0, or dBp/d^ 0, or hoth, these conditions can he readily 
evaluated, since i|r is known at the boundaries. 




4 . _ 


» 


■» 
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Analysis Problem 

The application of the characteristic equations to the analysis 
problem is quite simple and straightforward. In this problem the 
boundary curves are known, and the flow conditions are usually speci- 
fied at the inlet to the axisymmetric channel. A step-by-step solution 
along characteristic lines is then made, following the method of com- 
putation previously discussed, until the velocities throughout the field 
or channel have been obtained. From interpolation along the characteris- 
tic lines, the velocities across the outlet may then be obtained. In 
lieu of specified flow conditions at the inlet, the flow conditions may 
be specified at the outlet and a solution obtained for the flow in the 
channel and across the inlet. If, in a particular problem being solved, 
d(ru)/d^ = 0 and dEp/dijr = 0, integration to obtain the stream function 

along a characteristic line can be omitted if the streamlines are not 
desired; even in this case, however, concurrent integration of equa- 
tion (24) as the problem is being computed is recogmiended, since com- 
parison with the known values of the stream function at the boundaries 
furnishes an excellent check. 


Design Problem 

In applying the characteristic equations to obtain a solution to 
one type of design problem, all inlet conditions and part of the outlet 
conditions are specified and the velocities are prescribed where it is 
possible along a given wall, in order to obtain a solution far the 
opposite wall and the channel velocities. In the previous discussion 
in which figure 2 was considered, it was implied that the boundary con- 
ditions are known on a curve C containing such points as P]_ and Pg 

and that the solution is then extended into a region which contains 
such points as Pj and in which ru(ijr) and Hip(i)f) are known. On 

curve C, since w, v, ru, prp, and Erp are known, all the other vari- 
ables, such as d(ru)/di|r, dHp/dijr, \jr, u, p, and a 2 , maybe computed. 

Curve C should correspond to the line AB in figure 4 on which the 
inlet flow conditions are specified or to the given boundary curve 
BCDE. If the inlet flow conditions are specified, however, the flow is 
completely determined in the region bounded by the inlet AB, the por- 
tion of the given bo undar y curve BC, and the characteristic curve 

and the velocities may not be prescribed along BC. Conversely, if it 
is desired to prescribe the velocities completely along the curve BCDE, 
then only one of the velocity components v or w may be prescribed 
along AB. For the prescribed outlet conditions, the same considera- 
tions are valid in the region DEF. If the values of (ru) and Hr 
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are specified functions of i|r, various combinations of the remaining 
variables can be specified along the boundaries for the design problem 
as shown in the following table (see fig. 4): 


Case 

Given boundary 

Inlet 

Outlet 


BC 

CD 

DE 

AB 

EF 

1 



v,v 

w or v 
or v/w j 

w or v 
or v/w 

2 




w,v 

w or v 
or v/w 

3 

w,v 


v/w 

w or v 
or v/w 


4 

VjW 


v/w 

V;V 

v,v 


In case the inlet boundary falls along a curve AB, where B is 
between C and D, the inlet section AB has the same type .of boundary 
conditions as the outer contour BD; that is, the conditions on AB 
and ED are independent of each other. Similar relations apply at the 
outlet for an outlet line EE, when E is between C and D. A more 
general discussion of the relation between boundary shapes and boundary 
conditions can be found in reference 3. 

After the boundary conditions have been specified, a step-by-step 
solution is made along the characteristic lines until the velocities 
throughout the channel have been determined. Integration of the veloci- 
ties in the channel along characteristic lines gives the stream func- 
tion Points of equal i|r corresponding to the- full mass flow yield 

the final streamline or the required opposite boundary. 


NUMERICAL EXAMPLES 

The characteristic equations were applied and solutions were 
obtained for three examples of supersonic flow with whirl through axi- 
symmetric channels. These examples embody the analysis problem for 
irrotational inlet flow and the design problem for both irrotational and 
rotational inlet flow. The purpose of these solutions is to present 
some typical problems in order to illustrate the application of the 
characteristic equations. For all examples, the inlet flow is isoener- 
getic, that is, dEp/dijr = 0. The dimensions r and z. are made dimen- 
sionless by the inlet radius of the outer contour in al 1 examples. Also, 
the velocities, the static pressures, and the densities are made dimen- 
sionless by the inlet stagnation conditions of the velocity of sound, the 
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pressure, and the density, respectively. The stream function \|r is 
made dimensionless hy the value of i|r for the full mass flow. Dimen- 
sionless quantities are indicated hy the superscript *. The value of 
the ratio of specific heats x was taken equal to 1.400 for all 
examples. 


Example I 

The first example was the solution of the analysis problem for the 
irrotational flow through a channel representative of a vaneless diffuser 
for a supersonic compressor. 

Prescribed conditions . - The outer wall of this channel wa6 speci- 
fied by the equation 

- 2.2857Z* 2 - 0.49725r* 2 + 8.6790z* + 20.906r* - 8.99301^* - 20.408 = 0 

(27) 

between r* = 1.0 and r*= 1.9351. The inner contour was specified 


from z* = 0 to z* = 0.59552 by 

z* 2 + r * 2 - 6 . 0600r* + 3.9259 =0 (28) 

and from z* = 0.59552 to r* = 1.9351 by 

z* 2 - 0.4721z* + r* 2 - 4.3050r* + 2.7746 = 0 (29) 

The dimensionless inlet conditions were 

Vortex constant, r*u* 0.8640 

Axial velocity, w* 1.3377 

Radial velocity, v* 0 


Results . - The network for solution, or the points in the flow field 
of example I at which a solution was obtained, is shown in figure 5(a). 
Convergence of the or the T]- lines in this figure indicate the for- 
mation of a weak shock. If the |- or q-lines intersect, a solution 
can not be obtained at these points since the solution is multivalued. 

The shocks occurring are considered to be weak and the assumption is 
made that conditions remain isentropie in this region. The dashed lines 
in figure 5(a) show the streamlines of the flow through the channel 
obtained from the values of ljr* found during the solution. Figure 5(b) 
shows the contours of the resultant Mach number M throughout the 
channel which were obtained from the solution. This figure shows that 
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the convergence of the characteristic lines running from r* = 1.07, • 

z* = 1.1 on the inner wall to the point r* = 1.45, z * = 1.23 on the 
outer wall is accompanied hy a slight compression wave, and the reflected (k 

compression wave of this system from the outer to the inner wall is 
preceded hy an expansion wave (probably the result of the wall curvature) 
which shows up on figure 5(a) as a divergence of the characteristic 
lines. Also to be noted are the high outlet Mach numbers (which vary 
from 2.4 to 2.75) as compared with the inlet values (which vary from 
2.3 to 2.9), despite the conversion of tangential kinetic energy into 
pressure. Figure 5(c) shows the contours of the meridional component 
of the Mach number throughout the channel; this particular component 
emphasizes the effect of the pressure waves, since these waves affect 
only the meridional velocity components. This figure brings out the 
particular compression wave cited in figure 5(a) and also the fact that 
the reflected compression is preceded by an expansion wave. 

The contours of the tangential component of the Mach number through 
the channel are shown in figure 5(d) . This figure has been included to 
complete the plot of Mach number components. 

In figure 6 is shown the geometric relation between the a-coordinate 
used in figure 7 and the r*- and z* -coordinates. The a-coordinate 
is used for plotting the solution along the boundaries so that a compari- 
son between the outer and inner contours can be made more readily. The 
distribution of the static pressure, the resultant velocity, the meridi- 
onal velocity component along the outer and inner contours , and the area 
variations throughout the channel are shown in figure 7. 

The variations of the resultant velocity and the pressure on the 
outer and inner boundaries correlate with each other as might be 
expected, but the average of the values on the two walls shows only a 
small over -all change from inlet to outlet. On the other hand, the 
meridional component shows a steady increase (except for an initial drop 
on the inner contour) in the region of constant area (fig. 7) except 
near the end, where the decreasing area tends to keep the average meridi- 
onal component more nearly constant. The explanation of this anomalous 
behavior lies in the one -dimensional equation for mass flow 

W = pAV m _ (30) 

Differentiating this equation for the case of constant area after ^m 
has been expressed in terms of V and u yields 


V dV - u du 





(31) 
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The differentials of the density p and the velocity V are expressed 
in terms of pressure "by the Berno ull i equation 

V dV = ^ (32) 

and the isentropic relation 

— = constant (33) 

P r 

The differential of u is expressed in terms of r by use of the irro- 
tationality condition 


ru = constant 


The result obtained is 


d In p 
d In r 


rM e 

1 


(34) 


(35) 


For increasing radius, diffusion will require that this derivative be 
positive; hence, must be less than 1. 


Since for example I % is always greater than 1.0, an over-all 
decrease in pressure is to be expected from equation (35). 

Another significant factor in figure 7 is the fluctuation of the 
resultant velocity V*. The drop in velocity on the outer wall for 
45° < a < 50° shows the effect of compression waves originating on the 
inner wall at 35° < a, < 40° and gives an effect not calculable from 
simple one-dimensional flow or streamline flow with flow curvature esti- 
mated from wall curvature . There are similar important fluctuations 
elsewhere. In fact, the whole pattern of the resultant velocity from 
40° to 67° on the inner wall appears to be a reflection of the outer 
wall pattern from 28° to 55°, and the trough on the outer wall from 47° 
to 60° is a reflection of . the same pattern on Hie inner wall from 35° 
to 50°. 


Example II 

The second example is a solution of the design problem for the 
irrotational flow through a channel which might be representative of a 
transition channel between the outlet of a supersonic compressor and a 
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vaned diffuser. Such a channel might he desirable from the standpoint 
of changing the velocity distribution leaving the outlet of a supersonic 
compressor to obtain a flow better adapted for inflow to the diffuser 
blades. This transition channel might also be used to accomplish part 
of the required diffusion by contracting it to the extent that the 
Kantrowitz contraction ratio is reached for the meridional component of 
the Mach number. (The Kantrowitz ratio is the ratio of inlet area to 
throat area when a normal shock at the inlet would cause choking at the 
throat . ) 

Prescribed conditions . - The outer contour of the channel for this 
example was specified constant at a radius of 1.0 (dimensionless). The 
inlet conditions were the same as those specified in example I. The 
velocities could not be prescribed along the outer contour without over- 
prescribing the problem until the point (r* = 1.0, z* = 0.46, fig. 8(a)) 
at' which the characteristic originating at the inlet hub intersected the 
casing. From this point (r* = 1.0, z* - 0.46) until the end of the 
channel, the velocity was prescribed by the following equation 

w* = 15.396Z* 3 - 26.812z* 2 + 14.905z* - 1.3437 (36) 

This equation corresponds to diffusion of the meridional Mach number 
down to the value of the Kantrowitz ratio for an infinitesimal stream 
tube along the outer wall at a passage length of z* = 0.70, The other 
specified condition is that the radial, component of velocity is zero at 
the outlet. 

Results . - The solution of example II is shown graphically on fig- 
ure 8. The propagation of effects along characteristics (fig. 8(a)) 
again indicates the inadequancy of the one-dimensional or streamline 
curvature methods of computation. Between ijr* = 0 and ljr* = 0.2, for 
example, a throat exists at z* ~ 0.44. A similar minimum is seen 
between other streamlines. These minimums are located on the character- 
istic line between the points z* = 0.375, r* = 0.76 to z* = 0.70, 
r*= 1.00. Convergence and divergence of the characteristics indicate 
compression and expansion, respectively, as may be seen by comparison 
of figures 8(a) and 8(c). 

The constant velocity on the outer contour from z* = 0 to 
z*- = 0.46- (fig. 8(e)) is a consequence of the dependence of the velocity 
there on the inlet conditions (see characteristic curves, fig. 8(a)). 

(The possibility of a short diffuser of this type is thus eliminated.) 

The prescribed velocity drop from this point on was attained by com- .. 
pression waves formed by the calculated curvature of the inner contour 
from z* = 0 to z* = 0.375 (fig. 8(a)). The result is a pressure 
rise on the inner contour in the same region (fig. 8(e)). The outlet 
conditions of v* = 0 demand a reflex in this curvature and a pressure 
drop in the following region. . .. .... . 
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Example III 

Prescribed conditions . - The third example is similar to the second 
with the following differences: 

(1) The inlet tangential velocity was prescribed constant ins tead 
of inversely proportional to the radius and equal to 0.8640. The f lew- 
is therefore rotational; whereas, for cases I and II the flow is irro- 
tational. The inlet velocity vector at the outer contour was the same 
as that of example H, and .likewise, the radial velocity was zero. 

(2) The meridional velocity at the outer contour was prescribed 
after the intersection point with the characteristic line originating 
at the inner contour and the inlet. The value of w* in this region 
was: 

w* = 22.1308Z* 3 - 39.4467Z* 2 + 22.6930z* - 2.9145 (37) 

The initial and final values for w* in this region are the same in 
examples II and III. 

Results . - Because the rotational velocity is lower near the inner 
contour than for example H, the axial velocity components are higher 
and most of the differences between the two examples are a result of this 
velocity difference. The results are summarized in figure 9. Because 
of the expected decrease in Mach angles shown in figure 9(a), this 
diffuser is inherently longer than that for example II, but the diffusion 
rate at the outer contour was increased to maintain the same over-all 
axial length. All the effects noted in example II may be seen in exam- 
ple III with the modification for smaller values of £ (characteristic 
angles). The shortening of the region of diffusion on the inner contour 
and the shifting of the velocity rise region to the left, for example, 
may thus be explained from the change in the characteristic directions. 


SUMMARY OF RESULTS 

The axisymmetric flow between two surfaces of revolution with vor- 
ticity and tangential velocity was set up for solution by the method of 
characteristics. The criterion for the existence of characteristics is 
that the meridional component of the velocity be greater than sonic. 
Although the flow configuration is apparently simple, a holograph solu- 
tion similar to that for plane flows is not available because of the 
nonhomogeneity of the equations and because the coefficients of the dif- 
ferential equation are functions of position as well as of velocity. 

Some of the main results of this work are contained in the following 
statements : 
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1. In the numerical example of flow between two highly curved sur- 
faces of revolution, pressure waves arose which impinged on the opposing 
wall and reflected. Pressure variations therefore have a period struc- 
ture and a similarity will exist on the two boundaries. 

2. A one- dim ensional analysis showed that in the case in which 

along the flow path the flow area is constant and the radius increases 
so that the tangential component of velocity decreases, there will be a 
diffusion process only if M^< 1.0, where is the meridional com- 

ponent of the local Mach number. This fact was confirmed by the two- 
dimensional analysis of the example in which there was a radial flow. 

3. In all examples computed, the pattern of , convergence or diver- 
gence of the characteristic curves indicated compression or expansion 
waves which were the same as would be expected in two-dimensional plane 
flow. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 14, 1952 



APPENDIX - SYMBOLS 


The following symbols are used in this report: 
area 

sonic velocity of gas 

function of z (see equation (25)) 

enthalpy 

unit vector parallel to curve for 0 = constant and z 
constant (i = Vr) 

unit vector parallel to curve for z = constant and r 
constant (j = rV0) 

unit vector parallel to curve for r = constant and 0 
constant (k = Vz) 

local Mach number 

points defined hy figures 2 and 3 

static pressure of gas 

normal distance from axis of symmetry 

rotational component of gas velocity 

resultant gas velocity (V = iv + ju + kw) 

radial component of gas velocity 

axial component of gas velocity 

abbreviations (see equation (16)) 

abbreviations (see equation (17)) 

distance measured parallel to axis of symmetry 
angular position coordinate (see fig. 6) 
ratio of specific heats 
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slope of curve for | variable and rj constant 
slope of curve for T] variable and \ constant 

Tj one of the characteristic coordinates 

9 meridional angle (for cylindrical coordinates) 

| one of the characteristic coordinates 

p gas density 

a either r) or f in general form of equat ions 

i)r stream function (* = prjXV) 

Subscripts: 


i inlet 

m meridional component 

r,z,cr partial differentiation with respect to variable indicated 

by subscript 

I stagnation conditions 

9 tangential component 

1,2,3, 4 point or line positions (see figs. 2 to 4) 

Superscript : 

* dimensionless quantity 


REFERENCES 

1. Vavra, M. H. : Steady Flow of Nonviscous Elastic Fluids in Axially 

Symmetric Channels. Jour. Aero. Sci., vol. 17, no. 3, March 1950, 
pp. 149-156. 

2. Courant, R., and Friedrichs, K. 0.: Supersonic Flow and Shock Waves. 

Interscience Pub. Inc ., (New York) , 1948. 

3. Goldstein, Arthur W. : Axisymmetric Supersonic Flow In Rotating 

Impellers. NACA IN 2388, 1951. 



NACA TO 2768 


21 


4. Brand, Louis: Vector and Tensor Analysis. John Wiley and Sons, Inc., 

1948, p. 190. 

5. Southwell, R. V.: Relaxation Methods in Theoretical Physics. Claren- 

don Press (Oxford), 1946. 

6. Emmons, Howard W. ; The Numerical Solution of Compressible Fluid Flow 

Problems . NACA TO 932, 1944. 

7. Stanitz, John D. ; Two-Dimensional Compressible Flow in Turbomachines 

With Conic Flow Surfaces. NACA Rep. 935, 1949. (Supersedes NACA 
TN 1744.) 

8. Wu, Chung -Hua; A General Through-Flow Theory of Fluid Flow With Sub- 

sonic or Supersonic Velocity in Turbomachines of Arbitrary Hub and 
Casing Shapes. NACA TO 2302, 1951. 



22 


NACA TN 2768 



Figure 1. - Curved stream surface of revolution showing coordinates and velocity 

c opponents. 
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Figure 3. - Sketch to illustrate discussion of relations between characteristic 
lines end boundaries for supersonic analysis problem in annular channel. 
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(a) Coti tours of constant £ f Tf, and f M . 
Figure 8. - Solution far flow In example II* 
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(b) Contouro of oonfltant M and f*. 

Figure fl. - Continued. Solution for flow In example II. 
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Axial length, k* 

(d) Contjouro of oonotant K s mid f*. 

Figure 8. - Continued, Solution for flow In example U, 
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